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. ; ABSTRACT 

We have observed four low-luminosity active galactic nuclei classified as Type 1 LIN- 
Q . ERs with the X-ray Telescope (XRT) and the Ultra Violet-Optical Telescope (UVOT) onboard 

i-{ • Swift, in an attempt to clarify the main powering mechanism of this class of nearby sources. 

C/3 ' Among our targets, we detect X-ray variability in NGC 3998 for the first time. The light 

curves of this object reveal variations of up to 30% amplitude in half a day, with no signifi- 
cant spectral variability on this time scale. We also observe a decrease of ^30% over 9 days, 
with significant spectral softening. Moreover, the X-ray flux is ^40% lower than observed 
^ ' in previous years. Variability is detected in M 81 as well, at levels comparable to those re- 

\Q ' ported previously: a flux increase in the hard X-rays (1-10 keV) of 30% in ~3 hours and 

, variations by up to a factor of 2 within a few years. This X-ray behaviour is similar to that 

0^ ■ of higher-luminosity, Seyfert-type, objects. Using previous high-angular-resolution imaging 

data from the Hubble Space Telescope (HST), we evaluate the diffuse UV emission due to the 
Q>^ . host galaxy and isolate the nuclear flux in our UVOT observations. All sources are detected 

' in the UV band, at levels similar to those of the previous observations with HST. The XRT 

Q^ I (0.2-lOkeV) spectra are well described by single power-laws and the UV-to-X-ray flux ratios 

■ are again consistent with those of Seyferts and radio-loud AGNs of higher luminosity. The 

similarity in X-ray variability and broad-band energy distributions suggests the presence of 
similar accretion and radiation processes in low- and high-luminosity AGNs. 

5_j ■ Key words: galaxies: active — galaxies: nuclei — ultraviolet: galaxies — X-rays: galaxies 

' Facility: Swift, Chandra 



1 INTRODUCTION 

Low-luminosity active galactic nuclei (AGNs) are a common phe- 
nomenon, with a large fraction of all massive galaxies displaying 
some weak activity that is likely of non-stellar origin (see, e.g., 
Ho 2008, and references therein). Understanding the demograph- 
ics and physics of these objects is a necessary step towards the 
comprehension of super massive black hole activity in the local 
Universe and its past evolution. Low-Ionization Nuclear Emission 
Line Regions (LINERs), in particular, are a class of low-luminosity 
AGNs defined on the bas is of their optical spectral line ratios 
( iHeckmarj 1 1980l : IhoI liooi ) and can be further divided into vari- 
ous subclasses according to their different properties (Chiaberge 
et al. 2005; Gonzalez-Martin et al. 2009). A fundamental question 
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about these sources is the origin of their optical spectrum and multi- 
wavelength emission in general: it is not clear which fraction of the 
powering source is non-stellar and, if it is due to accretion, what 
is the regime of the accretion and the efficiency of the radiation 
conversion. It has been argued that the lack of X-ray variability 
(e.g., Roberts et al. 1999; Komossa et al. 1999; Georgantopoulos 
et al. 2002; Pellegrini et al. 2003; Ho 2008) , the non detecti on of 
broad Fe Kq lines down to stringent limits jPtak et alj|2004l) . and 
the weakness or absence of the characteristic 'big blue bump' in 
their optical/ne ar-ultra-vio let (UV) spectra, traditionally observed 
in Seyferts ( Ouataert et al.f l999; Chiaberge al. 2006), indicate that 
their engines may be intrinsically different from those of the more 
luminous AGNs, and could consist of radiatively inefficient accre- 
tion flows (RIAFs). 

However, X-ray observations of most of these sources have 
been sparse and not sensitive or not long enough to detect sig- 
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nificant variability, ex cept in M 81, one of the closest and best 
studied bona fide (e.g. iKewlev et ai]|2006l) LINERs. Furthermore, 
UV monitoring with the Hubble Space Telescope (HST) Advanced 
Camera for Surveys (ACS) of a sample of 17 LI NERs has reveale d 
the presence of bright and variable UV nuclei dMaoz et"aill2005l) . 
By coupling these UV measurements with non simultane ous X- 
ray m easurements with AiSCA, Chandra, and XMM-Newton. lMw^ 
( 12007 ) has shown that the UV-to-X-ray flux ratios in LINERs are 
similar to those of much-more-luminous Seyferts. Thus, contrary 
to the common paradigm, rather than being qualitatively different 
from Seyferts, LINERs may instead be 'scaled-down' analogues 
of Seyferts, with similar emission mechanisms operating in both 
classes. 

While the arguments of iMaoj | |2007|) are supported by accu- 
rate photometry, the UV and X-ray data of his sample were not si- 
multaneous, with measurements often separated by years. LINERs 
exhibit variability, which, albeit not of very large amplit ude (by fac- 
tors o f up to a few in the UV on time scales of years: iMaoz et al.l 
l2005l) . may undermine the results derived from non-simultaneous 
observations. In this paper, we approach this problem by obtain- 
ing, for the first time, simultaneous X-ray and UV data, along with 
X-ray variability information, for a small sample of LINERs. The 
NASA satellite Swift is well suited to this task because it is a flex- 
ible and efficient facility for long simultaneous and accurate UV 
and X-ray monitoring. Results of this work have been presented in 
preliminary form in Romano et al. (2009). 

Throughout this paper the uncertainties are given at 90% con- 
fidence levels for one interesting parameter (i.e., Ax^ ~ 2.71), 
unless otherwise stated. The spectral indices are parameterized as 
Fi, oc v~", where (erg cm~^ s~^ Hz~^) is the flux density as a 
function of frequency u; we also use F = a + 1 as the photon index, 
N{E) oc E^^ (ph cm~^ s~^ keV~^). This paper is organized as 
follows. In Section 2 we describe our sample, observations, anal- 
ysis, and results obtained with the two Swift instruments; in Sec- 
tion 3 we construct the spectral energy distributions (SED); and in 
Section 4 we discuss our findings. 



2 SAMPLE, OBSERVATIONS, DATA ANALYSIS AND 
RESULTS 

Our sample consists of the brightest LINERs in the iMaoj ^20071) 
UV sample, but excluding M 87, because of its prominent jet, 
which would dominate the UV emission at the Swiyi* resolution. The 
sample includes M 81, NGC 3998, NGC 4203 and NGC 4579 (see 
Table [T]l, which are all type-1 LINERs (i.e., wit h detected broa d 
Ha components) according to the classification bv lHoetalJ ( [l997h 
These four objects form a rather homogeneous sample (see lMao3 
l2007h : all have similar UV and X-ray luminosities (lO"" to 10*^ 
erg s~^), radio loudness parameters (~ 100), optical-to-X-ray in- 
dices aox (~ 1), central black hole masses (10'^ to 10** * M0, three 
of which are measured directly from stellar or gas kinematics) and 
Eddington ratios (—4.5 to —5.5). 

The Swift archivj^ already contained data for two of these 
sources, M 81 and NGC 4203, which had been observed as 
Swift fill-in targets, and which we retrieved. M 81 was observed 
for 3 3.7 ks and has an X-ray Telescope (XRT, iBurrows et al.l 
l2005h spectrum of excellent quality. NGC 4203 had a good XRT 
spectrum, although somewhat under-exposed compared to M 81 



(5.3 ks). We observed NGC 3998 and NGC 4579 with Swift for the 
first time as Targets of Opportunity (ToO) for 27.4 and 20.8 ks, re- 
spectively. 

2.1 XRT 

Table |2] reports the log of the Swift/XKT observations used for 
this work. The XRT data were processed with standard proce- 
dures (XRTPIPELINE vO. 11.6), filtering and screening criteria by 
using FTOOLS in the Heasoft package (v.6.4). Given the low 
count rate of the sources during the respective observing campaigns 
(< 0.5 counts s^^), we only considered photon counting data (PC), 
and further selected XRT grades 0-12. No pile-up correction was 
required. The source events were extracted in circular regions cen- 
tered on the source, with radii depending on the source intensity 
(10-20 pixels, 1 pixel ~ 2'.'37), while background events were ex- 
tracted in source-free annular or circular regions, depending on the 
field. 

Spectra were extracted for each XRT observation, as well as 
for the cumulative observing campaigns. Ancillary response files 
were generated with XRTMKARF, and account for different extrac- 
tion regions, vignetting and Point-Spread Function (PSF) correc- 
tions. We used the vOlO spectral redistribution matrices available 
in the Calibration Database maintained by HEASARC. All spectra 
were rebinned with a minimum of 20 counts per energy bin to allow 
fitting within XSPEC (vl 1.3.2). 
We extracted XRT light curves from the same regions as for 
the spectra for source and background in the standard bands, 0.2- 
lOkeV (total), 0.2-1 keV (soft, S) and 1-10 keV (hard, H). For our 
analysis, we considered several time bins ranging from 120 s to the 
typical orbit duration (< 5800 s), and the minimum time bin was 
chosen to ensure that each light curve point would have at least 30 
source-plus-background counts. The exposure requirements were 
that the bins be at least 50 % exposed. In each case, the light curves 
were corrected for PSF losses, i.e., losses due to the extraction re- 
gion geometry, and bad/hot pixels and columns falling within this 
region. 

The host galaxy contribution in the XRT extraction region 
was evaluated using archival Chandra images of these LINERs ob- 
tained with the longest exposures with ACIS-S. We evaluated the 
counts in the following regions: i) Rl, an annulus with outer ra- 
dius equal to the XRT extraction radius for each object and inner 
radius 4", i.e., the radius including 99% of the Chandra PSF; ii) 
R2, a circular region of 4" radius; Hi) R3, an annular region of radii 
4"and 4'.'5. We assumed that Rl only contained host galaxy pho- 
tons, R2 only active nucleus photons, while R3 was used to evaluate 
the galaxy counts in R2 assuming a constant host galaxy intensity 
profile. The measured counts in R3 were scaled to the area of R2 
and added to the counts in Rl, thus obtaining an estimate of the host 
galaxy counts in the whole XRT extraction region. The host galaxy 
contribution to the X-ray emission of the four LINERS ranges from 
6% (NGC 3998) to 16% (NGC 4203), and was therefore subse- 
quently ignored. 
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Figure 1. SwiftfKRT background-subtracted light curves (count rate in counts s ^) and hardness ratios (i.e., 1-10 keV to 0.2-1 keV flux ratios) of M 81. The 
binning time interval corresponds to the orbit duration. Panels (a) to (e) report the 0.2-1 keV (top), 1-lOkeV (middle) light curves and hardness ratio curves 
(bottom) of each pointing (see Table 2). Start times (t = 0) correspond to: (a) 2005 Apr 21.033 UT, (b) 2005 Aug 25.056 UT, (c) 2006 Jun 24.004 UT, (d) 
2006 Nov 18.045 UT, (e) 2007 Feb 6.906 UT. Panel (f) reports the 0.2-1 keV (fiUed circles) and 1-10 keV (open circles) Ught curves of M 81 between 2005 
and 2007. Each point is the average of the flux measured during each pointing in that given band and each curve is normalized to its average (0.16 and 0.29 
counts s~i for the 0.2-1 keV and 1-10 keV curves, respectively). The time origin (t = 0) corresponds to 2005 Apr 21.0 UT. 
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Figure 2. Swift/XKT background-subtracted light curves (count rate in counts s"'^) and hardness ratios (i.e., 1-10 keV to 0.2-1 keV flux ratios) of NGC 3998. 
The binning time interval corresponds to the orbit duration. Panels (a) and (h) report the 0.2-1 keV (top), 1-lOkeV (middle) light curves and hardness ratio 
curves (bottom) of each pointing (see Table 2). Start times (t = 0) correspond to: (a) 2007 Apr 20.628 UT, (h) 2007 Apr 29.11 UT. Panel (c) reports the 
0.2-1 keV (filled circles) and 1-10 keV (open circles) light curves of NGC 3998 in 2007. Each point is the average of the flux measured during each pointing 
in that given band and each curve is normahzed to its average (0.11 and 0.22 counts for the 0.2-1 keV and 1-10 keV curves, respectively). The time origin 
(t = 0) corresponds to 2007 Apr 20.0 UT. 
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Figure 3. Swift/XRT background-subtracted light curves (count rate in 
counts s~^) and hardness ratios (i.e., 1-10 keV to 0.2-1 keV flux ratios) of 
NGC 4203: 0.2-1 keV (top), 1-lOkeV (middle) light curves and hardness 
ratio curves (bottom). The binning time interval corresponds to the orbit 
duration. The time origin (t = 0) corresponds to 2005 Dec 25.004 UT. 



M 81 and NGC 3998 display a few very rapid (minute time 
scale) flux changes of 30 to 60%. Although this exceeds the typ- 
ical level of XRT photometric stability, which is about 10% (as 
established on the basis of the flux stability of the standard source 
PSR 0540-69, Cusumano, priv. comm.), we note that these occur at 
low flux levels, so that we cannot assess their authenticity. There- 
fore, we have considered only fluxes integrated over time scales no 
shorter than one orbit for further variability analysis. 

Since the Swift satellite must re-settle on the target at every or- 
bit, the target's location within the XRT field of view may change 
from orbit to orbit, and occasionally the source flux is greatly un- 
derestimated because of the vicinity of a detector dead column that 
affects the PSF to various degrees. We have systematically checked 
our observations for this effect and have corrected for it (see Grupe 
et al. 2007). However, because of the low flux level of our sources, 
in many cases the PSF reconstruction is not fully satisfactory and 
the corresponding flux points have been excluded from our vari- 
ability analysis. 

The orbit-averaged light curves of M 81 and NGC 3998 in- 
dicate variability. In Figure Ic, Id and 2b we show the flux and 
hardness-ratio time series of M 8 1 and NGC 3998 during the point- 
ings when significant variations were observed. The M 8 1 hard X- 
ray light curves in June (Fig. Ic, middle panel) and November (Fig. 
Id, middle panel) 2006 have ^ 30 (for 5 de grees of freedom, 
d.o.f.) and — 41 (for 19 d.o.f.) with respect to the average flux, 
respectively, corresponding to probabilities of constancy of 10~^ 
or less. The variations in hard X-rays have a maximum amplitude 
of 30% in a time scale of 3 to 12 hours. The variability in soft X- 
rays is less pronounced, but well correlated with the hard X-rays, 
so that the hardness ratio is not significantly variable on these time 
scales. The soft and hard X-ray light curves of NGC 3998 on 29 
April 2007 (Fig. 2b) are well correlated and have — 35 (13 



Figure 4. Swift/XKT background-subtracted light curves (count rate in 
counts s"'^) and hardness ratios (i.e., 1-10 keV to 0.2-1 keV flux ratios) 
of NGC 4579 during the 4 Swift pointings at this source: 0.2-1 keV (top), 
1-10 keV (middle) light curves and hardness ratio curves (bottom). The bin- 
ning time interval corresponds to the orbit duration. The time origin (t = 0) 
corresponds to 2007 May 15.039 UT. 



d.o.f.) and x — 44 (13 d.o.f.) with respect to the average flux, 
respectively, corresponding to constancy probabilities of 10~^ or 
less. The largest variability amplitude of NGC 3998 is about 30% 
in both X-ray energy ranges, on a time scale of about 12 hours. Nei- 
ther NGC 4203 (Fig. 3) nor NGC 4579 (Fig. 4) vary significantly 
in flux or hardness ratio on the orbit or longer time scale. 

We have evaluated the variability of M 8 1 and NGC 3998 dur- 
ing the individual Swift pointings also through the excess variance 
parameter af.^^, defined as in Nandra et al. (1997) and Turner et 
al. (1999; see also Ptak et al. 1998). This never exceeds the value 
of 0.02 for both sources and for both energy ranges, and therefore 
it is almost one order of magnitude smaller than would be expected 
- based on the X-ray luminosities of these two LINERs - from the 
inverse proportionality relation between X-ray luminosity and vari- 
ability established for Seyfert galaxies by Nandra et al. (1997) and 
Turner et al. (1999)0 

For M 81 and NGC 3998, which were observed on more than 
one epoch, with intervals larger than 1 day, we have evaluated the 
average flux at each epoch. The associated errors are the root mean 
squares of the averages, i.e. the standard deviations of the data 
points within the epoch, divided by the square root of their num- 
ber. The corresponding light curves are shown in Figures If and 2c 
for M 81 and NGC 3998, respectively. The long term variations are 
conspicuous. The historical flux of M 81 decreases between 2005 
April and 2006 June almost achromatically and then increases by 
2007 February by a factor of 2 in the soft X-rays and slightly less 



Note that the durations of our XRT observations are similar to the moni- 
toring times adopted by those authors, making the comparison of the (t^^^ 
parameters meaningful. 
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in hard X-rays. The variability on year time scales is well corre- 
lated in the two bands, with no significant variations of the hard- 
ness ratio. The 1-10 keV flux of NGC 3998 decreases by 30% in 
the nine days separating our two observations, while the 0.2-1 keV 
flux decreases more slowly, implying a softening of the spectrum. 
Indeed, the hardness ratio (i.e., the 1-10 keV to 0.2-1 keV flux ra- 
tio) decreases significantly between the 2 epochs, varying from HR 
= 2.088±0.004 to HR = 1.75±0.04. However, a spectral fit with a 
single power-law of each observation of this object does not show 
variability in the spectral parameters within the 90% confidence 
level. 

The mean XRT spectrum for each object was fit with a sim- 
ple absorbed power-law model with free absorption and photon 
index. Figure 5 shows the mean spectrum for each object in the 
sample, along with its best-fit model, while the fit results are re- 
ported in Table [T] The spectra show no significant absorption fea- 
tures superimposed on the power-law continuum; however, M 81 
shows hints of an Fe Ka emission line (6.4 keV, fixed), with an 
equivalent width EW = 640^410 eV. The F-test probability for 
this feature is 8.515 x 10""^, which corresponds to a 2.6-cr de- 
tection. NGC 4579 also shows a marginal (F-test probability of 
3.207 X IQ-^ or 2.1-cr) Fe Kq line at 6.4l:° j keV, with EW 
= 5501310° eV. The fluxes of these lines are (lOlf) x 10"" and 
(5^2^) X 10"^* ergcm"^ s~^, respectively. The equivalent widths 
are, withi n their large errors, broa dly consistent with the valu es re- 
ported bv lPellegrini et al. ! ('2000a') and Terashima et alj j2002l) . 

Although the uncertainties on the fitted hydrogen column den- 
sities (Table [T] column 4) are not small, all sources - except 
NGC 4203, whose spectrum has a poor signal-to-noise ratio - show 
significant evidence for absorption in excess of the Galactic one 
(also reported in Table[Tll. 

2.2 UVOT 

The S wift Ultra Violet-Optical Telescope (UVOT, iRoming et al.l 
l2005l) observed the four targets with the filters u (3465 A) uvwl 
(2600 A), uvm2 (2246 A), uvw2 (1928 A) simultaneously with the 
XRT. The filter choice was driven by the objective of maximizing 
the nuclear signal, which is dominant in the UV, while minimizing 
the stellar emission from the bulge populations of the target galax- 
ies. For each object, after verifying that the UVOT counts show no 
significant variability, we have coadded the images. 

The data analysis was performed using the UVOTSOURCE task 
included in the latest Heasoft software. This task normally cal- 
culates the magnitude by means of aperture photometry within a 
circular region of 5"-radius, which includes more than 99 % of 
the signal of a point source, and applies specific corrections due 
to the detector characteristics. Since the adoption of this standard 
aperture often resulted in over-subtraction of the host-galaxy back- 
ground, we used instead an aperture of 2" radius, and applied the 
proper correction factor between a 2"-radius aperture to a 5"-radius 
aperture (see below). Because of its proximity (3.6 Mpc), the host 
galaxy of M 81 occupies most of the field-of-view of UVOT, mak- 
ing it difficult to estimate the sky background to be subtracted from 
the nucleus and the galaxy light. The background was therefore 
computed from a 20'.'7-radius circular region separated by 6.' 5 from 
the nucleus. For uniformity, the remaining three LINERs, albeit 
more distant, were treated in the same way. 

Using archival HST images of these LINERs at 2500 A and 
3300 A ob tained with the H igh Resolution Channel (HRC) of the 
HST ACS jMaoz et al.l200^ . we have evaluated the nuclear-point- 
source-subtracted host galaxy contribution for each object in a 2 "- 



radius circular area centred on the nucleus, and have interpolated or 
extrapolated the UV host galaxy fluxes to the central wavelengths 
of the UVOT ti-band and UV filters, in order to subtract them from 
the observed fluxes. 

While the aperture-summed galaxy light is clearly detected in 
the HST images, the surface brightness is too low and the images 
are too oversampled to determine the surface brightness profiles 
of the galaxies on the scales of the apertures we use here. Different 
galaxy profiles will require different corrections for PSF losses out- 
side the aperture. For the galaxy light subtraction, we have there- 
fore considered two extreme cases: one in which the galaxy light is 
highly centrally concentrated, and one in which it has constant sur- 
face brightness. In the former scenario, the adoption of a 2 "-radius 
for the UVOT photometry results in a galaxy light loss similar to 
that suffered by the nuclear light, and therefore we first corrected 
the observed 2 "-radius flux to a 5 "-radius aperture, and then sub- 
tracted the galaxy flux evaluated from the HST images. In the latter 
case, where the light profile is flat, the galaxy light suffers no net 
PSF losses in the UVOT photometry, and we therefore first sub- 
tracted the HST galaxy flux from the total observed UVOT flux 
and then applied the 5"-radius aperture correction to the remain- 
der, which represents the point-source nuclear flux. 

The host galaxy contribution is dominant in the u band (it 
is generally comparable to the total flux observed by UVOT) and 
decreases with decreasing wavelength to the uvw2 filter (1928 A) 
where it contributes, at most, 50 % of the observed flux. The first 
galaxy-subtraction method described above predictably results in 
systematically larger nuclear fluxes than obtained with the second 
method. As also expected, the differences between the results of the 
two methods decrease as one goes to shorter wavelengths, where 
the galaxy light becomes less dominant. 

All UVOT magnitudes have been converted into fluxes us- 
ing the latest in -flight flux calibration factors and zero-points 
jPoole et al.ll2008 ). The average observed fluxes with their statis- 
tical errors are reported in Table [3] where we also list the host 
galaxy fluxes measured in the HST images and converted to the 
wavelengths of the UVOT filters, and the final nuclear fluxes ob- 
tained with the two subtraction methods described above. 



3 SED CONSTRUCTION 

We combined the XRT spectra, corrected for the total hydrogen 
absorption (Table [T] col. 4), and the UVOT fluxes, corrected for 
their host galaxy backgro unds and dere ddened using the A_b extinc- 
tion values compiled bv lMaozlJioOTh from Schlegel, Finkbeiner, 
& Davis (1998), and adopting the Galactic extinction curve of 
Cardelli, Clayton and Mathis (1989). These UV-to-X-ray spectral 
energy distributions (SED) are shown in Fig. 6. The errors asso- 
ciated with the UVOT nuclear fluxes are the sum in quadrature 
of the uncertainties associated with the observed UVOT fluxes 
[including the statistic al error r eported in Table [3] and the sys- 
tematic errors given in IPoole e t al. (2008)] and the uncertainties 
of the galaxy HST measurements (Table [3}, increased by 5 % - 
summed in quadrature - to take into account the uncertainty in 
the background subtraction on the HST images and the interpola- 
tion/extrapolation errors. Whenever the galaxy flux evaluated from 
the HST images exceeded the total flux observed by UVOT, we 
have used the total uncertainty to compute a 3-cr upper limit, which 
is indicated in the figure. 

For comparison, we have plotted also t he prev i ously measured 
UV-to-X-ray measurements [see Fig. 1 in iMaozl ( l2007h l and the 
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Figure 5. Swift/XKI spectra of the sample. The stepped curves represent the single absorbed power-laws which best fit the spectra (see Table[T]for spectral 
parameters). 



mean SEDs of radio-loud and radio-quiet quasars from lElvis et al.l 
( Il994t) . normalized to go through the geometric mean of the two 
nuclear fluxes at 1928 A obtained with the two galaxy subtraction 
methods described above. Note that the fluxes in this band suffer the 
least galaxy contamination, and therefore represent a solid lower 
limit to the intensity of the UV component. 



4 DISCUSSION 

We have observed four bright LlNERs with XRT and UVOT on- 
board Swift, to study their variability and their simultaneous UV-to- 
X-ray ratios, in an attempt to probe their emission mechanisms. In 
particular, high-luminosity AGNs ~ Seyferts and QSOs - are gen- 
erally thought to be powered by geometrically thin and optically 
thick accretion disks. Our new data can bring new insight to the 
question of whether low-luminosity AGNs - specifically LlNERs 
- are simply lower-power analogues of Seyferts and QSOs, or are 
instead powered by some fundamentally distinct process, such as a 
RIAE 

Among our results is the first detection of X-ray variability 
in NGC 3998, in a range of time scales from a few hours to days. 



with a ~30% amplitude. This is in contrast to earlier reports of 
no X-ray variability during Bep poSAX and XMM -Newton obser- 
vations (Pellegrini etal. 2000b: Pt^ et alj|2004l) . The reason for 
this difference may be related to the larger sensitivity of XRT 
with respect to BeppoSAX, and to XRT's higher scheduling flex- 
ibility, which allows for longer moni t oring p eriods than possible 
with XMM-Newton. IPellegrini et alj ( l2000 tf) have cited the ab- 
sence of variability in NGC 3998 as evidence that the accretion 
is advection-dominated in this source. Other workers (Ptak et al. 
1998; Ho 2008) have, in general, invoked non-variability in X-rays 
as an indicator of a RIAF mode. Our finding of significant short- 
term X-ray flux variations in a low-luminosity LINER points, in- 
stead, to a similarity with whatever accretion mode is occurring 
in higher luininosity AGNs. However, the excess variance parame- 
ter of the NGC 3998 XRT light curves does not obey the inverse 
relation between X-ray luminosity and variability amplitude de- 
termined for high-luminosity AGNs (Nandra et al. 1997; Turner 
et al. 1999) when this is extrapolated to the X-ray luminosity of 
NGC 3998 (1.1 X 10^" erg s"^ in 2-10 keV). Instead, the vari- 
ability amplitude is similar to that typical of Seyferts of luminosity 



10* 



10 erg s (2-10 keV). The range of validity of the re- 
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Figure 6. Spectral energy distributions of our sample based on the simultaneous SwiftrXKT and UVOT data from the present campaign (filled symbols) and 
from previous, non-simultaneous, HSTW and Chandra, XMM- Newton o r ASC A X-ray o bser vations [empty squares; Ho et al. ( 2001); Terashima et al. ami); 
IPtak et al.l <2004l) ; lLa Parola et al.l <2004l) ; lMaoz et alj <2005l) ; ICap pi et al. (2006); Maoz (2007)]. The optical-UV data, coiTected for Galactic reddening, result 
from the observed UVOT fluxes, minus the HST host galaxy fluxes assuming that the host galaxy profile is strongly centrally concentrated (filled triangles) or 
constant with radius within the UVOT PSF (filled circles). When the latter method yields a negative flux (see Table 3), we replace both nuclear fluxes with a 
3-0" upper limit. The X-ray data are corrected for both Galactic and intrinsic absorptio n. The power laws fitting the XRT (solid bow-tie) and archival (dotted 
bow-tie) X-ray spectra are reported with their 90% confidence boundaries. Simil arlv tolMaoj j2007h , we have superimposed on the data the mean spectra of 
radio-loud (long-dashed curve) and radio-quiet (dot-dashed curve) quasars from lElvis etalJ h994), and we have normalized them to the geometric mean of 
the 1928 A nuclear fluxes obtained with the two methods of galaxy subtraction. We also show lines that connect the 2500 A and 2 keV fluxes (solid) and the 
lines predicted by the S06 relation between ciox and the monochromatic UV luminosities (short dash). 
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lation may thus not extend to the low X-ray luminosities of the 
objects studied here. 

A second LINER in our sample, M 81, has already been shown 
to be X-ray variable on many time scales, from hours to years (Pel- 
legrini et al. 2000a; La Parola et al. 2004; Page et al. 2004; Markoff 
et al. 2008; lyomoto & Makishima 2001; Young et al. 2007). Our 
findings confirm this behaviour and set this object, along with 
NGC 3998, as another X-ray-variable low-luminosity AGN that is 
similar in this respect to high-luminosity AGNs. It has been sug- 
gested that the mechanisms active in this obj ect are similar to thos e 
that cause the hard states in X-ray binaries if Markoff et al.ll2008h . 
Our observation of spectral softening accompanying source bright- 
ening on a time scale of years is indeed reminiscent of the behavior 
exhibited, on a shorter time scale, by Galactic sources. 

Our XRT timing results on NGC 4203 and NGC 4579 are in- 
conclusive. NGC 4203 exhibits no significant X-ray variability on 
any time scale, despite being, among the sources in Maoz ( 2007j), 
the one that varied with the largest amplitude in the UV. One rea- 
son may be that our monitoring period was limited, and hence we 
did not catch the stochastically occurring variations. Alternatively, 
for this object the errors associated with orbit-averaged points are 
somewhat larger than those of M 81 and NGC 3998, so that inter- 
orbit variations of amplitude similar to the inter-orbit variations in 
M 81 and NGC 3998 would not be detected at high significance. 
A prominent, variable Fe Ka emission line has been detected by 
ASCA and XMM-Newton in NGC 4579, and has given rise to var- 
ious interpretations as to its origin from a standard accretion disk 
or a truncat ed disk toewangan et al.ll2004l ; iTerashima et al. I ll998l. 
I200d l2002l). Our weak XRT detection of the feature does not add 
significantly to this debate. 

Our sources are all well detected in at least one UV filter 
(Fig. 6). T he UV- t o-X-r ay SEDs compare rather well with those 
reported bv lMaoj ( l2007h and obtained with non simultaneous UV 
and X-ray data. The X-ray spectrum of the LINERs always lies 
about a factor of 3 to 10 above the normalized radio-loud quasar 
template, depending on the object and on the wavelength. The 
UV flux level has not changed dramatically in any of the sources 
(Fig. 6), alth o ugh th ese are among the most UV-variable sources 
in the iMaozl j2007h sample. We note, however, that the UVOT 
fluxes are more uncertain than the HST measurements, owing to 
the galaxy subtraction and to the lower angular resolution and pho- 
tometric accuracy of UVOT. The XRT fluxes and spectra are also 
similar to the older ones, except for NGC 3998, which XRT has 
detected in a state lower by ^^40% than observed in 1999 and 2001 
by BeppoSAX and XMM-Newton, respectively (Ptak et al. 2004, 
see also Fig. 6b). 

We have evaluated the Eddington ratios of our LINERs us- 
ing the bolometric lu minosit i es an d the central black hole mass 
estimates compiled in iMaozl ( |2007|) . The bolometric luminosities 
were obtained fro m the normalized templates of radio-loud AGNs 
( lElvis et al.|[l994l) from 3465 A to 0.2 keV and the observed X-ray 
spectrum from 0.2 to lOkeV. In fact, the template may be a more 
reliable description of the spectrum at UV wavelengths (consider- 
ing that some of our UV points are upper limits), while at X-rays 
the template always underestimates the real flux (see Fig. 6). In 
the radio-loud AGN spectrum, the UV-to-X-ray flux dominates the 
SED, representing more than half of the total emission. Moreover, 
considering that the intrinsic AGN spectra probably lack the con- 
spicuous observed bump cent red at 1 /xm, due to r eprocessed radi- 
ation at shorter wavelengths ( iMarconi et al. I l2004l) . our computed 
bolometric luminosities may adequately represent the intrinsic lu- 
minosities after all. The Eddington ratios, reported in Table [T] are 



about 0.5 dex higher than those computed bv iMaozl ( l2007t) . who 
only included the UV luminosity, but are still low, as typical for 
LINERs. 

By normalizing the AGN templates to the UV emission, we 
note that the precisely determined XRT spectral slopes follow the 
average radio-loud AGN shape, although, as noted, the spectral nor- 
malizations exceed the prediction by 0.5 to 1 order of magnitude 
(Fig. 6). Strateva et al. (2005) and Steffen et al. (2006, hereafter 
S06) have studied the dependence of the optical-to-X-ray colour 
index, ctox, of AGNs (defined - using our notation for the spec- 
tral indices, see last paragraph of Section 1 - as aox = -0.3838 
log(/2feeV /I ° ), where lohev and / • are the rest-frame X- 

cv^fcev/ 2500a" ^'^''^ 2500A 

ray and UV monochromatic luminosities in erg s^^ Hz~^, respec- 
tively) on the monochromatic rest-frame UV luminosity, and have 
found an empirical linear relation between those quantities (the pa- 
rameters of the relation in the two works coincide, within the er- 
rors; we have used that of S06, which is based on a larger sam- 
ple). We have extrapolated Uox in the S06 relation, Uox = 0.137 
log (/ ") - 2.638 (note that we have changed the sign of their 

■= ^ 2500 a' ^ ^ " 

Eq. 2, to be coherent with our spectral index notation), by two or- 
ders of magnitude to low UV luminosities, corresponding to the 
observed monochromatic UV luminosities of the LINERs in our 
sample (extinction-corrected, see Table 1, col. 9). This extrapola- 
tion is plotted as a short-dashed line in Fig. 6. It can be compared 
to a power law passing through our ohsen'ed luminosities at 2500 
A and 2 keV, shown in Fig. 6 as a solid line, with errors omitted, 
for clarity. 

In NGC 3998 and NGC 4579 the observed values of aox co- 
incide, to within the errors, with the extrapolation to low luminosi- 
ties of the S06 relation. In M 81 and NGC 4203, the less UV- 
luminous sources in our sample (see Table 1), aox.ots has an in- 
termediate value, similar to those of Seyferts, and the extrapolated 
Uox is significantly flatter than the observed one and predicts X-ray 
luminosities that are 4-5 times higher than observed. As already 
noted by S06 based on their data alone, and further pointed out by 
Maoz (2007) based on his non simultaneous UV and X-ray mea- 
surements, the S06 relation probably flattens at luminosities below 
~ 10^^ erg s~^ Hz~^. Overall, our present simultaneous UV and 
X-ray photometry strengthens the case for a limiting value of Uox 
below some critical luminosity typical of Seyferts. 

In summary, for the LINERs under study, the significant ul- 
traviolet emission, the X-ray intraday variability, and the X-ray 
and multiwavelength spectral similarity with radio-loud AGNs all 
point to continuity and similarity with AGNs of higher luminos- 
ity. In view of this, it is not clear that distinct accretion and radia- 
tion mechanisms are required in the different luminosity regimes. 
Further Swift/XKI and UVOT monitoring, or more optimally, new 
HST high angular resolution UV observations of larger samples, 
accompanied with XRT simultaneous coverage, could shed more 
light on these questions. 

We note, finally, that we cannot exclude that, in these sources, 
the emission at X-ray and/or UV bands is partially due to a weak 
jet, some evidence of which has been reported for all of our four 
sources, either directly (for M 81, based on resolved radio struc- 
ture detection and radio polarization) or indirectly (for NGC 3998, 
NGC 4203, NGC 4579, based on a flat radio spectrum, a com- 
pact and variable radio core, and high brightness temperature) from 
high resolution radio imaging (Bietenholz et al. 2000; Bower et 
al. 2002; Filho et al. 2002; Anderson et al. 2004; Bietenholz et 
al. 2004; Ros & Perez Torres 2008). In the case of NGC 3998, 
a jet has been proposed by Ptak et al. (2004) as possibly respon- 
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sible for the multiwavelength emission. Jetted high energy radi- 
ation in this LINER would be compatible both with the rapid 
flux variations and with the frequency-dependent variability am- 
plitude that we observe ov er 9 days (Fig. 2c) , which is typical of 
non-thermal sources (e.g. lUlrich et al.lll997l) , and which is dis- 
similar to the X-ray variability observed in Seyferts [usually of 
larger amplit u de at s ofter energies, e.g. , Matsuoka et al. (199C); 
Nandra et al. ( 1990); Edelson et al. ( 2000); Turner et all fcOOll) ; 
Uttlev & M cHardv (200 5); Terashima et al. (2009)1 and in weak- 
line radio-galaxies re.g. jGliozzi et alT l 20081) 1. However, it would 
be premature to draw a conclusion on the exact nature of the X-ray 
emitting process from the detection of spectral softening accom- 
panying flux decrease between only two epochs in this individual 
LINER. More work on the correlated radio and X-ray variability in 
LINERs - in both the observational and the theoretical directions - 
is necessary to assess the role of a jet in the multiwavelength emis- 
sion (see e.g. Brenneman et al. 2009). 
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Table 1. Target parameters. 



Name 


Distance 
(Mpc) 


(1020 cm-2) 


(1020 cm-2^ 






pi 

(0.2-1 keV) 


(1-lOkeV) 






M81 


3.6 


5.55 


10.45«:«3 


2.04 ± 0.04 


0.966/290 


11.2 ±0.1 


14.9 ±0.4 


1.4 


8.8 X lO"*' 


NGC 3998 
NGC 4203 


13.1 
15.1 


1.01 
1.11 


6.85t;;2g 

^- '*-2.74 


1.95 ±0.06 

-I O-1+0.24 


0.948/188 
0.699/16 


6 2+0-1 
"•^-0.8 

-1 4+0-2 


9 8+0'5 

^■"-0.4 


6.0 
4.0 


1.1 X 10^5 
1.1 X 10^'' 


NGC 4579 


21 


2.97 


7 (,,+1.43 


1.92 ±0.07 


0.972/145 


4.1 ±0.1 


6.8 ±0.4 


9.8 


7.4 X 10-5 



Hydrogen column densities derived from lKalberla et al] j2005l) and consistent with the Ag extinction values reported in jMaoj2007h . using a typical Milky 
Way gas-to-dust ratio, 5 X 102i cm" 2 mag" 
^ Hydrogen column densities from the XRT spectral fits. 

Photon index, oc . 

Unabsorbed flux in units of 10~^^ erg cm~2 §—1 
" Unabsorbed monochromatic luminosity at 2500 A, in units of 1025 erg s^ ^ Hz^ ^ . 
' Eddington ratio based on the observed SEDs. 

Table 2. XRT observation log. 



Name 


Sequence^ 


Start time (UT) 


End time (UT) 


Exposure 






(yyyy-mm-dd hh:mm:ss) 


(yyyy-mm-dd hh:mm:ss) 


(s) 


M81 


00035059001 


2005-04-21 00:47:46 


2005-04-21 23:11:33 


1563 




00035059002 


2005-08-25 01:20:04 


2005-08-25 23:59:59 


5022 




00035059003 


2006-06-24 00:05:06 


2006-06-24 08:02:57 


4140 




00035059004 


2006-11-18 01:05:11 


2006-11-20 23:54:55 


20986 




00259527001 


2007-02-06 21:44:17 


2007-02-06 23:13:04 


1973 


NGC 3998 


00030916001 


2007-04-20 15:03:58 


2007-04-20 16:49:57 


1198 




00030916002 


2007-04-29 02:39:06 


2007-04-30 14:32:27 


26216 


NGC 4203 


00035477001 


2005-12-25 00:05:03 


2005-12-25 07:04:59 


5261 


NGC 4579 


00030939001 


2007-05-15 00:56:39 


2007-05-15 14:07:01 


7712 




00030939002 


2007-05-16 01:08:52 


2007-05-16 23:39:56 


7870 




00030939003 


2007-05-17 01:18:44 


2007-05-17 12:30:58 


3490 




00030939004 


2007-05-18 10:53:03 


2007-05-19 23:59:23 


1743 



^ We only considered data collected in PC observing mode. 
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Table 3. SwiftlUWOT Observations. 



Name 


Band 


Wavelength(A) 


Observed flux''''' 


Galaxy flux"'"^ 


Nuclear flux(l)'*'<l'* 


Nuclear flux(2)''-'«'* 


M 81 




3465 




^.J-U u.uo 


1.2 ± 0.6 


—6.3 ± 0.6 






2600 




-L.UU ^Z U.UO 


3.4 ± 0.4 


1.5 ± 0.4 




null 


2246 




u.oo HI u.uo 


1.8 ± 0.4 


0.9 ± 0.4 




uvw2 


1928 


3.39 ±0.29 


0.26 ±0.03 


3.1 ±0.5 


2.7 ±0.5 


NGC 3998 


u 


3465 


5.73 ±0.02 


3.9 ±0.1 


1.8 ±0.7 


1.0 ±0.7 




uvwl 


2600 


3.13 ±0.01 


2.4 ±0.1 


0.8 ±0.3 


0.3 ±0.3 




uvm2 


2246 


2.46 ±0.01 


1.82 ± 0.03 


0.6 ± 0.4 


0.3 ± 0.4 




uvw2 


1928 


3.11 ±0.01 


1.39 ± 0.03 


1.7 ±0.4 


1.4 ± 0.4 


NGC 4203 


u 


3465 


2.44 ± 0.07 


2.35 ± 0.03 


0.1 ±0.3 


-0.4 ± 0.3 




uvwl 


2600 


1.42 ±0.04 


0.56 ±0.03 


0.9 ±0.1 


0.8 ±0.1 




uvm2 


2246 


0.62 ±0.02 


0.27 ± 0.03 


0.3 ±0.1 


0.3 ±0.1 




uvw2 


1928 


0.96 ± 0.04 


0.13 ±0.03 


0.8 ±0.1 


0.8 ±0.1 


NGC 4579 


u 


3465 


3.11 ±0.06 


2.74 ±0.03 


0.4 ±0.3 


-0.2 ±0.3 




uvwl 


2600 


1.84 ± 0.03 


1.33 ±0.03 


0.5 ±0.1 


0.3 ±0.1 




uvm2 


2246 


1.55 ±0.06 


0.92 ± 0.03 


0.6 ±0.2 


0.4 ± 0.2 




uvw2 


1928 


1.51 ±0.03 


0.62 ± 0.03 


0.9 ±0.2 


0.8 ±0.2 



Not corrected for Galactic reddening, in 10 erg s 'cm '. 

Total flux, evaluated in a 2 "-radius circular area centred on the nucleus, and scaled to the equivalent flux in a 5 "-radius area using the UVOT PSF. The 
errors are statistical only (1-a). 

Extranuclear galaxy flux within a 2 "-radius circular area centred on the nucleus, obtained by interpolating/extrapolating from archival HST ACS HRC 
images at 2500 and 3300 A. The errors are statistical only (l-cr). 
Net nuclear flux, obtained by assuming that the host galaxy radial profile is highly centrally concentrated. 

Net nuclear flux, obtained by assuming that the host galaxy radial profile is nearly flat, i.e. constant within the UVOT 2 "-radius aperture. 
The errors are propagated from the Swift UVOT and HST ACS HRC uncertainties and include the systematic uncertainties. 



